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Abstract. Unbonded Braced Frame (UBF) system is one of the efficient damping 
systems. It has been made of a smart combination of steel and concrete or mortar. In 
this system, steel bears axial force of earthquake-induced as compressive or tension 
forces without strength decreasing. Concrete or mortar around the steel core acted with 
decreasing of brace slenderness as a constraint for brace and prevents of brace 
buckling during seismic axial load as compressive or tension forces. In this paper, a 
steel structure with UBF system was modeled using Finite Element Method (FEM) 
which Modal and Nonlinear Time-history Analyses (NTA) were utilized considering 
the effects of Near-field and Far-field earthquakes. Furthermore, three various Optimal 
Sensor Placement (OSP) methods were used and Genetic algorithm was selected to act 
as the solution of the optimization formulation in the selection of the best sensor 
placement according to structural dynamic response of the UBF system. Results show 
that with a proper OSP algorithm for Structural Health Monitoring (SHM) in the UBF 
structures can be detected weak and sensitive points in comparison without utilizing 
UBF system. 

Keywords: Unbonded Braced Frame (UBF), Finite Element Method (FEM), Optimal 
Sensor Placement (OSP), Nonlinear Time -history Analysis (NTA), Genetic Algorithm 
(GA) 

1. Introduction 

Providing reliable mechanisms for dissipation of the destructive earthquake energy is a key 
for the safety of structures against intense earthquakes. Inelastic deformations can limit the 
forces in members allowing reasonable design dimensions; and provide hysteretic energy 
dissipation to the system. 

An Unbonded Brace is a buckling-restrained brace where Euler buckling of the 
central steel core is prevented by encasing it over its length in a steel tube filled with mortar. 
The basic principle in the construction of the most popular type of unbonded brace is to 
prevent Euler buckling of a central steel core by encasing it over its length in a steel tube 
filled with concrete or mortar. The term "unbonded brace" derives from the need to provide a 
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slip surface or unbonding layer between the steel core and the surrounding concrete, so that 
axial loads are taken only by the steel core. The materials and geometry in this slip layer must 
be carefully designed and constructed to allow relative movement between the steel element 
and the concrete due to shearing and Poisson's effect, while simultaneously inhibiting local 
buckling of the steel as it yields in compression. The concrete and steel tube encasement 
provides sufficient flexural strength and stiffness to prevent global buckling of the brace, 
allowing the core to undergo fully-reversed axial yield cycles without loss of stiffness or 
strength. The concrete and steel tube also helps to resist local buckling [1]. 

Most Structural Health Monitoring (SHM) research to date has focused either on 
global damage assessment techniques using structural dynamic responses, or on limited local 
independent damage detection mechanisms. The idea of using dynamics characteristic data 
for structural parameter identification and damage detection is especially attractive because it 
allows for a global evaluation of the condition of a structure. 

The performance of the system varies depending on several factors including sensor 
or sampling locations. Location of sensors is one of the most important factors in a 
monitoring network and it needs to be optimized to maximize system performance and 
reduce the cost of the system [2]. Hemez and Farhat [3] extended the effective independence 
method in an algorithm where sensor placement was achieved in terms of the strain energy 
contribution of the structure. Miller [4] computed a Gaussian quadrature formula using the 
functional gain as a weight function, and thought that the nodes of the quadrature formula 
gave the optimal locations for sensors. Hiramoto et al. [5] used the explicit solution of the 
algebraic Riccati equation to determine the optimal sensor/actuator placement for active 
vibration control. Wouwer et al. [6] presented an optimality criterion for the selection of 
optimal sensor locations; the criterion was based on a measure of independence of the sensor 
responses. Worden and Burrows [7] used a number of different methods to determine an 
optimal sensor distribution based on the curvature data. Also, Cobb and Liebst [8] and Shi et 
al. [9] have reported the Optimal Sensor Placement (OSP) for the purpose of detecting 
structural damage. In recent years, Cruz et al. [10] and Yi et al. [11] performed new 
researchers about OSP for two different structures. In this research, a steel structure with 
UBF system was modeled using Finite Element Method (FEM) which employed Modal 
Analysis (MA) and Nonlinear Time -history Analysis (NT A) considering the effects of Near- 
field and Far-field earthquakes. Furthermore, three various OSP methods were utilized and 
Genetic Algorithm (GA) was selected to act as the solution of the optimization formulation in 
the selection of the best sensor placement according to structural dynamic response of the 
steel building. 

2. Structural Modeling Description 

In this paper, a 6-story structure was chosen for three-dimensional analysis with a regular 
plan. The structure under the study consists of Moment Resisting Frame (MRF) in one 
direction and uses UBF system in another direction. The braced frame and moment frame 
were modeled in parallel connected by a rigid diaphragm. A 3D view of the FE structural 
model is shown in Figure 1 . 

This building considered has the same plan and the LRFD design method was utilized 
in the designing procedure. It was located in the high damage risk zone of the Tehran city. 
Beam-to-column connections were considered rigidly and ST-37 as steel material property 
was used for steel elements. Table 1 presents the material properties, in which E, d, G are 
modulus of elasticity, Poisson ratio and shear modulus, respectively. 




Figure 1 3D view of the FE structural model 
Table 1 The material properties data 



Yield strength 
(f y -kg/cm 2 ) 



Ultimate tensile 

strength 

(f„-kg/cm 2 ) 



E (steel- kg/cm 2 ) 



i9 (steel) 



G (steel- kg/cm 2 ) 



2400 



3600 



2.1 x 10° 



0.3 



7.84193 x 10 



3. Analysis Procedure 

3.1 Nonlinear Time-History Analysis (NTA) 

Six sets of recorded earthquake ground motions including near-field and far-field seismic 
records were selected from the Pacific Earthquake Engineering Research database [12] for 
the purpose of dynamic analysis. Nonlinear time-history analyses were conducted using a FE 
software to evaluate dynamic response of presented structural system. Non-linear dynamic 
analyses using 6 scaled earthquake ground motions were performed for each model. Seismic 
responses are obtained in x, y and z directions. The time histories obtained are maximum 
base shear, moment and lateral displacements in the x and y directions. It should be noted that 
due to model accounts directly for effects of material inelastic response, the resulting internal 
forces are reasonable approximations of those expected during the design earthquake. The 
evidences in Tables 2 and 3 have been attained from nonlinear dynamic analyses of the 
presented brace system in the 6-story building considering near-field and far-field 
earthquakes. 



Table 2 Maximum results of nonlinear time-history analyses for UBF system in the 6-story building 

using near-field earthquakes 
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Design base 


Design base 


Top-story 


Top-story 


Earthquake 


shear (Kgf) 


shear (Kgf) 


moment (Kgf 


moment (Kgf 


displacement 


displacement 




in the x- 


in the in- 


-m) in the x- 


-m) in the y- 


in x-direction 


in y-direction 




direction 


direction 


direction 


direction 


(cm) 


(cm) 


Chi-Chi 


833650.55 


588334.4 


4009351.56 


5493207.08 


8.81 


9.26 


(1999) 














Chi-Chi 


711204.14 


720344.24 


8793468.76 


6442753.11 


8.38 


8.94 


(1999) 














Northridge 


208266.52 


283474.42 


1436272.9 


1386389.85 


8.22 


8.51 


(1994) 















Table 3 Maximum results of nonlinear time-history analyses for UBF system in the 6-story building 

using far-field earthquakes 

Design base Design base Design base Design base Top-story Top-story 

Earthquake shear (Kgf) shear (Kgf) moment (Kgf moment (Kgf displacement displacement 
in the x- in the y- -m) in the x- -m) in the y- in x-direction in y-direction 



direction 



direction 



direction 



direction 



(cm) 



(cm) 



Chi-Chi 
(1999) 


796083.58 


642755.62 


6355504.9 


8317546.93 


7.63 


7.92 


Chi-Chi 
(1999) 


832035.34 


588766.66 


4011572.56 


5498486.99 


7.28 


7.71 


Northridge 
(1994) 


408549.52 


146321.25 


1869778.25 


2735895.98 


7.14 


7.48 



3.2 Modal Analysis (MA) 

Natural time periods or natural frequencies are the important characteristics of a structure. It 
can be used to analyze the results obtained by dynamic analysis. To evaluate natural 
frequencies, modal analysis of UBF braced building has been carried out. Table 4 shows the 
first four natural time periods and maximum modal responses of the proposed FE model 
respectively. Figure 2 illustrates also the mode shapes of braced structure. 

Table 4 Modal results: periods and maximum displacement in the X, Y and Z directions 



Mode no. 


Period (sec) 


Max Dis. 
(X-cm) 


Max Dis. 
(Y-cm) 


Max Dis. 
(Z-cm) 


1 


1.12 


12.42 


11.46 


11.95 


2 


1.08 


9.44 


10.15 


10.17 


3 


0.74 


7.28 


8.44 


8.22 


4 


0.41 


5.32 


6.11 


6.06 





2" a Mode 





3 rd Mode 



4 m Mode 



Figure 2 The effective mode shapes of UBF steel structure 



4. Genetic Algorithm (GA) 

To find the optimal solution of the objective function given in the previous section, genetic 
algorithms were used. Genetic algorithm is well known due to its simplicity and robustness in 
the solution of complex problems, and these characteristics fit well to the problem described 
here [13]. The length of an individual chromosome depends on the number of sensors and all 
integer numbers in a chromosome should be unique. To place the sensors simultaneously, a 
genetic algorithm is an effective method. 

GAs have been widely used in sensor placement type problems [14]. They have been 
used to search for the optimal locations of actuators in active vibration control [15]. Another 
type of problem where GAs have been successfully utilized is the placement of vibration 
isolators to reduce the transmissibility of undesirable vibrations to an optical laboratory table 
[16]. However, these methods often produce some invalid strings in the evolution process. 

GA has been proved to be a powerful tool to OSP, but it also has some faults that 
need to be improved. For example, one sensor location may be placed where two or more 
sensors or sensor numbers are not equal to a certain number [17]. A good surveillance system 
must detect the aftermath of a contaminant spill as soon as possible with maximum 
reliability. Earlier detection gives more time for the system managers to react. Thus an 
objective of the optimization is to maximize the detection time that can be defined in an 
average sense and an objective function. 

In this optimization method, information about a problem, such as variable 
parameters, is coded into a genetic string known as a chromosome (individual). Each of these 
chromosomes has an associated fitness value, which is usually determined by the objective 
function to be maximized or minimized. Each chromosome contains substrings known as 
genes, which contribute in different ways to the fitness of the chromosome. It should be noted 
that an essential characteristic of a GA is the coding of the variables that describe the 
problem. A binary coding method can be used. The method is to transform the variables to a 
binary string of specific length. Also, Crossover is the operator that produces new individuals 
(offspring) by exchanging some bits of a couple of randomly selected individuals (parents). 
Furthermore, Mutation operates on a single individual with a small probability. With this 
operation, one or more bits are chosen at random from the individual and changed into a 
different symbol [18]. 

5. Sensor Placement Optimization Procedure 

Numerous techniques have been advanced for the OSP problem and are widely reported in 
this literature. These have been developed using a number of approaches and criteria, some 
based on intuitive placement or heuristic approaches, others employing systematic 
optimization methods. The sensor placement optimization can be generalized as "given a set 
of n candidate locations, find the subset of m locations, where m <n, which may provide the 
best possible performance" [11]. 

In the case under investigation the fitness function is a weighting function that 
measures the quality and the performance of a specific sensor location design. This function 
is maximized by the GA system in the process of evolutionary optimization. As known, the 
measured mode shape vectors in the SHM have to be as possible linearly independent, which 
is a basic requirement to distinguish, measured or identified modes. Hence, three different 
OSP algorithms were utilized which include Modal Assurance Criterion (MAC), Extended 
MAC (EM AC) and Transformed Time-history to Frequency Domain (TTFD). 



5.1 MAC Algorithm 

The simple way to check linear dependence of mode shapes is to calculate the MAC [19]. It 
is equivalent to maximize the angles formed by unit mode shape vectors, or to minimize the 
dot product between them, which is the same as the MAC. The MAC without mass weighting 
is just to compare the direction of two vectors. When two vectors lie in the same direction or 
near, the MAC value or the correlation coefficient is one or approximate one. Small 
maximum off-diagonal term indicates less correlation between corresponding mode shape 
vectors, and renders the mode shapes discriminable from each other. The fitness function 
presented in this section is constructed by the MAC, the biggest value in all the off-diagonal 
elements in the MAC matrix. The reason for the selection of this kind of fitness function is 
that the MAC matrix will be diagonal for an optimal sensor placement strategy so the size of 
the off-diagonal elements can be taken as an indication of the fitness. The MAC can be 
defined as Equation (1), which measures the correlation between mode shapes: 

MACu = — f- 1 — Kf — (1) 

where, <2>j and &j represent the ith andy'th column vectors in matrix <2>, respectively, and the 
superscript T denotes the transpose of the vector. In this formulation, the values of the MAC 
range between and 1 , where zero indicates that there is little or no correlation between the 
off-diagonal element MAC t j (i =£ j) and one means that there is a high degree of similarity 
between the modal vectors [19]. 

Then the MAC fitness function is given as Equation (2): 

/ = 1 - max (abs(MAC tJ )) (i * ;') (2) 

For this attempt, the size of the searching space is the number of nodes on the FE 
model excluding the constrained nodes and the vibration nodes of the selected modes. The 
MAC algorithm achieves this objective as follows. First, an intuition sensor set (much less 
than the required number of sensors) is selected based on experience and requirements of 
structural topology for visualization of mode shapes. Second, it adds other available 
candidate sensors one by one, and selects one that minimizes the maximum off-diagonal 
element of the MAC matrix at each step. Third, the MAC repeats the second step by adding 
one sensor at a time until a required number of sensors are selected. When performing the 
OSP method via GA technique, certain parameters are required such as; constraints, Fitness 
function, population, generation, crossover method, mutation rate etc. 

The locations selected for the 3 accelerometers needed to be installed. It should be 
noticed that due to the nature of the GA method, the results are usually dependent on the 
randomly generated initial conditions, which means the algorithm may converge to a 
different result in the parameter space. These values are very close to the optimal value. 

5.2 EMAC Algorithm 

To overcome the contra-decreasing problem of the original MAC algorithm, a forward- 
backward combinational extension is developed as follows by Dongsheng [20]. An EMAC 
algorithm proposed to overcome the disadvantages of traditional MAC algorithm with the 
introduction of a forward-backward combinational approach. First, an intuition sensor set, U 



(including, to say, a number of sensors, s ) is chosen. Then, one sensor is added to this initial 
set until a preset number of sensors, which is somewhat larger than the number of sensors as 
required, for instance, ten percent more than required (l.ls ), is reached. This is the same 
forward sequential MAC procedure. The extension differs from the original forward approach 
in the stopping criterion. The EMAC algorithm is continued to obtain a sensor set, [/ l5 
consisting a certain number of sensors (to say, s l5 s x < l.ls ) larger than the required one s 
where the original MAC stops [20]. 

Secondly, one sensor at each step is excluded from the sensor set f/ x until the required 
number of sensors s is reached. This is the backward sequential MAC approach, the 
essential extension to the forward one. Therefore, two function curves are established. One is 
the curve of the maximum off-diagonal term with respect to the number of sensors increasing 
from s to s-l be obtained in the first stage, and the other is the curve of the maximum off- 
diagonal term with respect to the number of sensors decreasing from s t to s found in the 
second stage. Both curves are compared and the one with a smaller value at the point s is 
selected. Which curve is to be selected, depends on the abilities of the forward and backward 
approaches to minimize the maximum off-diagonal terms of the MAC matrix. In this manner, 
the maximum off-diagonal term of the MAC matrix may, in many instances, be further 
minimized than the traditional MAC algorithm [20]. 

Naturally, the forward stopping number of sensors s x in the first step can be varied 
according to the structure under consideration. The effects of various numbers s t of sensors 
on the selection set (including 5 sensors) of the above two step processes can be compared 
and the one with the smallest maximum off-diagonal term of the MAC matrix can be chosen. 
This can be implemented as the third step, if necessary. 

As before section, 3 accelerometers needed to be installed. Figure 3, illustrates which 
the best fitness values tend to a constant quicker along as the number of generations increases 
despite many fluctuations occurring that are caused by the genetic operators of crossover and 
mutation's search procedures. 

There is one note about the influence of the choice of the intuition sensor set, U , on 
the final selection of sensor positions. If a newly added sensor conflicts with one or several of 
the original intuition set, the intuition set maybe reformed if the exclusion of certain sensor 
from the original intuition set U is not considered to be much detrimental to the mode shape 
visualization. Afterwards, the two steps can be recomputed. The optimal locations for the 
MAC and EMAC algorithms are illustrated in Table 5. 
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Figure 3 Results of optimization process using EMAC algorithm ("A" and "F" parameters are 
corresponding to Average and Fitness respectively). 



Table 5 Comparison of the optimal sensor locations of the MAC and EMAC algorithms 



Sensor No. 


1 


2 


3 


MAC algorithm 


1 


2 


5 


EMAC algorithm 


2 


4 


6 



5.3 TTFD Algorithm 

Due to lack of exact dynamic response on the previous presented OSP algorithms; in this 
research, authors introduced TTFD algorithm using time -history analysis. In this method, the 
results of NT A procedure are employed so that the optimal locations can be acquired exactly. 
At the first, the mentioned results in the time-domain should be transformed to frequency- 
domain. Any signal can completely be described in time or in frequency domain. As both 
representations are equivalent, it is possible to transform them to each other. This is done by 
the so-called Fourier Transformation (Equation (3)) and the Inverse Fourier Transformation 
(Equation (4)), respectively: 



U(ja>) 



= j u 



(t). e~ }(0 - c dt 



(3) 



u{t) 



-t-uo 

= j- [ U(jo)- e M c da> 



(4) 



After domain transformation, the periods obtained corresponding to Displacement-Time 
record of NTA. In following the minimum and maximum periods of the modal analysis 
(Table 4) were compared with the obtained periods by NTA procedure and so the maximum 
NTA displacements were corresponded to maximum modal responses. This procedure 
likewise was undertaken for the 2 n and 3 r modes. Figures 4 and 5 illustrate the Fourier 
Amplitude-Frequency diagram for far- field and near-field earthquakes respectively. 
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Figure 4(a) Response Acceleration-Time 
diagram of Chi-Chi Earthquake (far-1999) 



Figure 4(b) Time-Frequency domain 
transformation of Chi-Chi Earthquake (1999) 
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Figure 5(a) Response Acceleration-Time 
diagram of Northridge Earthquake (1994) 



Figure 5(b) Time-Frequency domain 
transformation of Northridge Earthquake (1994) 



The Complete Quadratic Combination (CQC) method was utilized as a well known 
modal combination technique in this paper but after Time -Frequency transformation authors 
replaced the equivalent values of NTA with the modal analysis results. Hence, all of the 
required parameters values in the CQC such as maximum modal responses in the nth and mth 
modes, rigid response coefficient etc were replaced with NTA transformed and equivalent 
results. According to Iranian seismic code [21] and also ASCE code [22], since maximum 
modal responses do not occur for different modes in an earthquake simultaneously therefore, 
it should estimate maximum modal responses in different members of structure via a 
statistical method. The mentioned statistical method should be as indicated by the maximum 
displacements of different modes and it should include the effects of probable interactions 
between different displacements close together that are derived from different modes. The 
peak value of the total response U is estimated by combining the peak modal response of 
individual modes using modified double sum equation; this is given by [23]: 



1\ iV — _L iV 



Pnm'-'n'-'m 



(5) 



n=l 



n=l m=n+l 



where, U n and U m are the maximum modal responses in the nth and mth modes, and p nm is 
the modified correlation factor defined as: 



= a n a m + [(1 - a n 2 ) (1 



i 

a m 2 )] 2 Pn 



(6) 



where, a n is the rigid response coefficient in the nth mode and p nm is the correlation 
coefficient of the damped periodic part of modal responses, given by the well known CQC 
rule. For damped periodic modes, a = 0, and modified double sum equation reduces to CQC 
and for p nm = 0, modified double sum method reduces this to the Square Root of Sum of 
Squares (SRSS). Equations (5) and (6) include the effect of rigid response of high frequency 
modes in the modified correlation coefficient p nm . The rigid response coefficient a n is 
defined as [23]: 



a„ = 



// X n {t)u g {t)dt 



t d a n x a 



Xrv u Q 



(7) 



where, X n {t) is the acceleration response, a n x and a u a are the standard deviations of X n (t) 
and u g (t) respectively and t d is the duration of responses. The modal responses, having a 



frequency less than rigid frequency also have a rigid content and the value of a gradually 
reduces from one to zero from a key frequency f 2 to another key frequency f t [23]. The key 
frequency f 2 is the lowest frequency at which the rigid response coefficient becomes 1 and 
the key frequency f t the highest frequency at which the rigid response coefficient becomes 
zero. An approximate equation for a n can be represented by a straight line between the two 
key frequencies f x and f 2 on a semi logarithmic graph, this is given by [23]: 

a n = -% < a n < 1 (8) 

where f n is the modal frequency (Hz) and the key frequencies f x and f 2 can be expressed as, 

r ->A max 

Jl = 2nS v max ( 9 ) 



f _ (/l + 2/ r ) 

h- * 



(10) 

3 

where, S A max = maximum spectral acceleration, S v max = maximum spectral velocity and f r 
= rigid frequency. 

Now the optimization procedure can be performed by GA. In this study, the number 
of variables is equal to 6 and also 18 constraints have been considered. The height of braced 
building has been divided into 6 ranges. It should be implied that the height of the structure is 
18.0 meters. So, the length of each range was found to be 3.0 meters. At the first, p nm was 
calculated then the square of total response U 2 was computed using p nm , U n and U m 
parameters. Each of the variations should be less than maximum modal response 
corresponding to that mode in the each range. Equation (4) was employed as Fitness function 
in this framework. This function was maximized by the GA system in the process of 
evolutionary optimization. When performing the OSP method via GA technique, certain 
parameters are required such as; constraints, Fitness function, population, generation, 
crossover method, mutation rate etc. Each of these runs started with a random initial 
population that was uniformly distributed within the same range. Figure 6, clearly shows that 
the best fitness values tend to a constant quicker along as the number of generations increases 
despite many fluctuations occurring that are caused by the genetic operators of crossover and 
mutation's search processes. Finally after obtaining U 2 values via optimization process, 
interpolation carried out in each range considering the initial results of the square of total 
response so acquired the optimal locations. 

According to this optimization method and considering GA results, the most proper 
locations for layout of smart sensors in the building have been presented as the final results in 
Tables 6 and 7. 
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Figure 6 Results of optimization process in all of stages, evolution of the best function value 

performance of optimized sensor layout as compared to stages of 1 to 6 ("A" and "F" parameters 

are corresponding to Average and Fitness respectively). 

Table 6 Details of optimization by GA in each 6 stages using a far-field earthquake 



Stage 


1 (0-3.0 m) 


2(3.01-6.0 
m) 


3 (6.01-9.0 
m) 


4(9.01-12.0 
m) 


5(12.01-15.0 
m) 


6(15.01-18.0 
m) 


X direction 
(m) 


2.44 


4.92 


8.32 


11.33 


14.03 


16.94 


Y direction 
(m) 


2.27 


4.32 


8.09 


11.12 


13.73 


16.38 


Z direction 
(m) 


2.08 


4.01 


7.96 


10.78 


13.02 


16.07 


Table 7 Details of opt 


imization by 


GA in each 6 


stages using a 


near-field earthquake 


Stage 


1 (0-3.0 m) 


2(3.01-6.0 
m) 


3 (6.01-9.0 
m) 


4(9.01-12.0 
m) 


5(12.01-15.0 
m) 


6(15.01-18.0 
m) 


X direction 
(m) 


2.56 


5.22 


8.61 


11.47 


14.23 


17.33 


Y direction 
(m) 


2.34 


5.03 


8.22 


11.21 


13.87 


16.94 


Z direction 
(m) 


2.68 


5.46 


8.73 


11.56 


14.52 


17.41 



6. Conclusions 



Researches show that application of health monitoring can detect weak points of structure 
and it can help prevent possible damages that may occur in the future, it also plays a key role 
on improving the structural performance. In this paper, three reliable methods for OSP of the 
UBF steel structure, based on modal and time-history analyses result, for the evaluation of 
objective function have been demonstrated. Real coded elitist genetic algorithms with 
uniform parent centric crossover operators and mutation operator have been developed and 
used for the implementation of the optimal placement. 

In this research, MAC, EMAC and TTFD methods were investigated for OSP 
procedure. MAC algorithm is a common method in this field whereas EMAC is a new 
algorithm that has been improved MAC which both of them utilizes the free vibration 
analysis results. Furthermore, a novel approach was proposed by authors for OSP which was 



adopted TTFD algorithm. This novel method uses time-history analysis results as an exact 
seismic response despite the MAC and EMAC algorithms which utilize modal analysis 
results. Therefore, for the health monitoring of a steel structure, layouts of smart sensors in 
the mentioned points will be used to reveal the structural condition with the appropriate 
layouts. These will be derived from numerical methods. Once the layout of smart sensors are 
placed in the mentioned areas of building a safe method for health monitoring of different 
structures can be implemented. 
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